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We construct a model of quantum metrology inspired by the computational model known as 
deterministic quantum computation with one quantum bit (DQC1). Using only one pure qubit 
together with l fully-mixed qubits we obtain measurement precision at the standard quantum limit, 
which is typically obtained using the same number of uncorrelated qubits in fully-pure states. The 
standard quantum limit can be exceeded using an additional qubit, which adds only a small amount 
of purity. We show that the discord in the final state vanishes only in the limit of attaining infinite 
precision for the parameter being estimated. 


Quantum information science has transformed how we 
view many computation, communication, and precision- 
measurement tasks. Emerging quantum technologies 
promise to solve problems that are intractable or im¬ 
possible using classical counterparts. However, in many 
cases the origins of quantum enhancements remain the 
subject of debate. Entanglement unambiguously plays a 
critical role in many tasks that use pure states, but this 
often ceases to be true when noise is added to the pic¬ 
ture [T. One of the most studied tasks that uses noisy 
qubits is provided by a model called DQC1, introduced 
by Knill and Laflamme [2j. DQC1 performs a specific 
type of classically-hard computation using highly-mixed 
quantum states, and thereby seriously challenges the no¬ 
tion that pure-state entanglement plays an essential role 
in quantum computation. 

The task performed by DQC1 is to estimate the nor¬ 
malised trace of a quantum circuit U that acts on a col¬ 
lection of l register qubits, as depicted in Fig. [lja). The 
initial state comprises one “clean” pure qubit together 
with register qubits that are maximally mixed, and only 
unitary gates are used for the computation. Remarkably, 
the precision of the estimate does not scale with the size 
of U. It is intuitively clear that DQC1 achieves an expo¬ 
nential speedup over any classical algorithm which finds 
and sums the 2 l eigenvalues for U, and there is now con¬ 
siderable evidence which supports the existence of a true 
quantum speedup for DQC1 [3]. Several works have also 
analysed how the computational power of DQC1 changes 
as resources, such as additional pure qubits and measure¬ 
ments, are added @HS], see Fig. 0b). 

Some studies have also investigated the role of en¬ 
tanglement and quantum discord [313 in the speedup 
achieved by DQC1 mm- It has been found that the 
discord generated at the output of DQC1 m , by uni¬ 
tary transformations which are randomly selected accord¬ 
ing to the Haar measure, remains a fixed proportion of 
the maximum possible as the unitary transformations in¬ 
crease in size. However, the amount of entanglement gen¬ 
erated by these unitary transformations is vanishing. It 
is not yet known that happens to entanglement or discord 
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FIG. 1. Circuits for DQC1. (a) A DQCl-complete problem is 
to compute the normalised trace of a unitary transformation. 
For the circuit, a Hadamard gate is applied to the control 
qubit followed by controlled-unitary transformation on the 
register, and measurement of the control. After several runs of 
the circuit, an estimate is obtained for ( a x )+i(a y ) = tr{U)/2 l . 
The precision of this estimate depends only on the number of 
runs of circuit. The protocol also works when the control 
qubit is partially pure at the start - as given by the state 
given in Eq. |5|. In this case, the number of runs must be 
increased by a factor 1/e 2 to achieve the same precision as 
when the control qubit is initially pure, (b) In the general 
DQC1 problem, an additional n ~ log(l) pure qubits can be 
introduced without altering the computational power of this 
model [4j|5 . 

at intermediate steps in a DQC1 computation. In con¬ 
trast, it well known for pure-state quantum computation 
that unbounded entanglement is necessary for exponen¬ 
tial speedups when using circuits composed of gates of 
bounded size [Ti. 

We now turn to quantum metrology, and the problem 
of achieving quantum advantage for precision in the task 
of phase estimation, which is used for highly-sensitive 
measurements of physical parameters 17211 11 1 . Phase- 
estimation strategies that cannot exploit quantum fea¬ 
tures are subject to the standard quantum limit (SQL) 
for precision, given by A<f> = l/\/n where n particles are 
used as probe and <fi is the parameter to be estimated. 
For example, this limit applies when n single pure qubits 
in the |+) = (|0) + | l ))/-\/2 state are used to measure the 
phase for a Pauli rotation 

M0 = e* 1 ^ 0 where g = |l)(l|. (1) 
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However, when a GHZ state |+ n ) = (|O) 0 ”-|-|l) 0n )/\/2 is 
used as the probe state with G = 1 0 ji the precision 

scales at the Heisenberg limit A <j> = 1 /n, which is the 
best precision achievable m- 

Inspired by DQC1, we ask whether a large ensemble 
of mixed qubits can be used as the basis of a powerful 
sensor. We consider a model where only one (or few) 
clean qubits are accessible, and only one qubit can be 
measured at the end |15) . Physical systems where our 
model is most relevant include NMR fT6l - fT8] and some 
cold-atom systems mm- For these systems often only 
bulk operations on the register qubits are available - 
which is so say the same operation is applied to every 
register qubit, optionally under global control. Hence, 
we add a bulk-operations constraint to our model. 

For mixed-state models of phase estimation, recent re¬ 
sults challenge any presupposed link between entangle¬ 
ment and quantum advantage for measurement preci¬ 
sion. Ref. m considers an algorithm for multi-parameter 
estimation using DQC1. This algorithm uses an adap¬ 
tive protocol based on a series of estimates with differ¬ 
ent interactions times, to achieve a final precision scaling 
with the inverse total interaction time. Ref. [22] analyses 
the situation where a unitary circuit is used to prepare 
probe states from n uncorrelated qubits in the state p e 
given below in Eq. ([5]). Strategies using probe states 
with only classical correlations [23], which is to say they 
are diagonal in the a z basis up to local-unitary transfor¬ 
mations |24j . are compared with strategies with exploit 
entanglement and quantum discord (defined as in [25]). 
It was found that circuits which generate non-classical 
correlations can achieve a quadratic quantum advantage 
compared to circuits generating only classical correla¬ 
tions at fixed e in Eq. 0). This result holds even for 
small values of e where there is no entanglement but large 
amounts of discord, and the amount of discord also grows 
with n. Another recent analysis considers phase estima¬ 
tion using an interferometer, where the spectrum of the 
interferometer Hamiltonian is fixed but not its eigenbasis. 
The authors found that the minimum amount of statisti¬ 
cal information that can be extracted about the unknown 
phase in the problem also constitutes a measure of the 
discord-type correlations in the probe state [25]. 

Parameter estimation .— In the classical theory for pa¬ 
rameter estimation m a probability distribution p is 
subjected to a process that is a function of a single pa¬ 
rameter cf>. The process alters the initial distribution p 
into p(()>), which depends on the value of </>. Differentiat¬ 
ing between the initial and the final distributions allows 
for the determination of the value of </>. The uncertainty 
in this value is bounded by the Fisher information, which 
is given by: 

A</> > —^=, with F = y ]2 (2) 
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FIG. 2. Illustration of our general scheme for phase estima¬ 
tion: There are n pure qubits, m semi-pure qubits, and l 
fully-mixed qubits (I 2 = 1/2). The first qubit is the control, 
and the remaining qubits constitute the register. Only bulk 
operations are permitted for gates used to prepare the probe 
state and implement the readout procedure. A bulk Cnot is 
used to prepare the probe, and then each qubit is subjected 
to the unitary operation given in Eq. 0 for which cj> is to be 
determined. The readout procedure is adaptive: 0(r) is the 
estimate for <j> after the first r— 1 rounds. This value is used 
to configure the readout circuit, which is a Cnot followed by 
controlled-v^r) on all qubits of the register, and measurement 
on the control. 


When using a quantum system, the initial and final 
probability distributions are replaced by density opera¬ 
tors cr and cr(</>) respectively. The final state is measured 
by a positive operator valued measure (POVM) {n^.} to 
yield classical probabilities pu = trplfccr^)], from which 
F in Eq. can be computed. If the process which is pa¬ 
rameterized by <j> is unitary, then the Fisher information 
when optimised over all POVMs is given by the quantum 
Fisher information [30] : 

^q = 4^ (A ;~y K^m-)i 2 , ( 3 ) 

. Aj t Aj 

i>j J 


where {A;} are the eigenvalues, {\4>i)} are the eigenvec¬ 
tors of a, and G is the Hamiltonian generator of the 
phase shift. This formula for Fq is very powerful: It 
yields the lower bound for the precision of </> without 
needing the explicit form of the optimising POVM, and 
enables a straightforward comparison between different 
initial states. 

The setup .— Our model uses three registers: one with 
n pure qubits; one with m qubits with finite purity as 
given in Eq. 0); and one with l fully-mixed qubits. 
Along with these three registers, there is one pure qubit 
in state |0) which serves as the control. The total initial 
state is 


with 


f 7o = |O)(O|®|O)(O| 0n 0pf TO 






0 < e < 1. 


(4) 

(5) 


and pk is the probability for observing outcome k. The (Later we will take the limit of n, m —> 0 to develop corn- 
above inequality is the Cramer-Rao bound [251 125] . pare our results with DQC1.) To prepare the probe state 
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a we apply the Hadamard gate to the control qubit fol¬ 
lowed by a Cnot gate for all qubits in the register. Next 
each qubit in the register is allowed to evolve freely un¬ 
der the unitary operation given in Eq. 0 . The readout 
procedure consists of another controlled operation and 
measurement of the control qubit. The full protocol is 
shown in Fig. [2j 

To compute Fq for a above, we note that <7o has eigen¬ 
vectors of the form |±; Bq ; B™; H[); here Bf represents 
a binary string of length a with Is appearing i times, and 
the semicolons separate the control qubit and the three 
registers. There are ("') ( l k ) such eigenvectors each with 
eigenvalue 

\ + = 2^7 + (6) 


when the control qubit is in state |+), and the eigen¬ 
value is XJ = 0 otherwise. After the first Cnot gate the 
eigenvectors are 


^ (|0; So”; Bf; B{) ± 11; C"; C$»; C{)) , (7) 


where C “ is the NOT of F“, i.e., | C 7 “) = erf 


B j)- 

IXIIx® 1 *. 


The generator of the phase shift is G = 
where is the identity operator on all but ;cth qubit, 
and x runs from 1 to n + m + l. Next, we note that the 
components of the eigenstate of the prepared state are 
eigenvectors G: 


G 10; BJ; B?*; B l k ) = (j + k ) |0; F 0 "; BJ-B l k ) , 
G |1; C£; Cf; C l k ) = (n + 1 + m - j + l - k) 


v 11 ■ r' n ■ f 171 • c 1 \ 

x I 1 ; o 0 , Cjj , O fc/ ) 


and therefore 


( 8 ) 

(9) 


Wp k '\GW,+ ) =-{j + k)(Bp\B™){Bi,\B l k ) 

± -(n + m + l — j — k + 1) 

x <c>|C3"><c£i|c£>. (io) 


Eq. (10) is non-zero only when j = j' and k = k'. We 
note that the numerator of the first term in Eq. (|3j) is 
the difference in two eigenvalues, and therefore it is only 
necessary to consider Hence, 




j =0 K J / k—0 

= l + m(l — e 2 ) + (1 + n + em) 2 . (11) 


We can make several observations concerning Eq. ( p~T| ) : 
(i) Fq is always greater or equal to the SQL value, which 
is 1+l+m+n. (ii) The SQL is attained when m = n = 0, 
i.e. the case which is analogous to DQC1 (31). (iii) If 
e is small (or even 0) there is a linear contribution of 
m corresponding to size of the register of partially-pure 


F 



FIG. 3. (Colour online.) Fisher information, F, computed 
from the probability distribution in Eq. (121 for n = 6 and 
m = l = 0 (dotted); n = l = 0 and m = 11 with e = 0.49 
(dashed); and l = 48 and n = m = 0 (solid). F is independent 
of cj) when register has no mixed qubits. All three cases have 
the same value for the quantum Fisher information Fq. 


qubits, (iv) (n + l) 2 exhibits the well-known quadratic 
enhancement for entangled pure states of n + 1 qubits, 
and there is an additional contribution equivalent to em 
extra pure qubits. 

Readout procedure. — Apart from preparation of the 
probe state, attention must be given to the bulk- 
operation requirements for implementing the measure¬ 
ments for the readout procedure. In other words, we 
need to consider how Fq given in Eq. © can be at¬ 
tained via a suitable POVM, which in general can re¬ 
quire entangled measurements [521 . For our model, the 
following method suffices (illustrated in Fig. [2]): a bulk 
Cnot gate is performed, followed by a bulk controlled- 
vg( r ) where vg( r ) = exp {— i9(r)a z }, and a measurement 
on the control qubit. 9(r) here is taken to be the es¬ 
timate of (f> following the (r — l)th round. The initial 
estimate 0(0) can assume no prior knowledge of <j>. In 
each successive round our estimate for is improved, 
i.e., \9(r)—<j>\ < \9(r— l)—(j>\, using an adaptive Bayesian 
update or maximum-likelihood method to maximize sen¬ 
sitivity [33] . 

The measurement of the control qubit along the a x 
direction yields probability distribution 

q ± (r) = - (1 ± x(r)) where (12) 

x(r) =Re{e i(n+lMr) cos z (u;(r)) 

x (cos(w(r)) + iesin(w(r))] m |, (13) 

and w(r) = 9(r) —<j>. The value for F computed from this 
probability distribution, using Eq. ([2]), yields a value that 
approaches the quantum Fisher information in Eq. 
as w(r) ss 0. That means that the adaptive protocol 
described above will yield the optimal Fisher information 
as the estimate 0(r) approaches </>. We have plotted three 
cases in Fig. [3] 

One+e clean qubit metrology. — It is clear that the 
metrology protocol presented in Fig. [2] is a special case 
of DQC1 computation, provided n + m ~ log(Z). In fact, 
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we can think of the circuit as an application of a bulk 
controlled-unitary operation where the unitary 

operation is 


/ \ 0(+m+rt 

U u (r) = \%ve(r)^ x u^,a x ) . (14) 

Using this feature we will examine the discord hypothesis 
for DQC1 due to [101 . 

However, let us first consider the case where n = 0, 
to = 1, and l > 0. In this case the probe state <r is 
entangled for any value of e > 0. This can be understood 
by noting that a non-positive partial transpose for a of 
the state results from applying a Cnot gate on the state 
in Eq. |5]), controlled on |+). Another way to see this is 
by noting that the value for Fq here beats the SQL TTT] : 

Fq = l + 2 + 2e > l + 2. (15) 


In other words, even with one qubit with finite purity we 
can attain better precision than what is possible classi¬ 
cally. Adding more qubits to the registers for initially 
partially-mixed and pure qubits, the entanglement (be¬ 
tween the control and registers) will increase as well as 
the value for Fq. 

One-pure-qubit metrology. — We now let n = to = 0, 
i.e., consider a l + 1 qubit state with only one pure qubit 
and l qubits in fully-mixed state. From Eq. © we see 
that Eq has the SQL value of l + 1 qubits, and that the 
SQL is attained using only one pure qubit and l fully- 
mixed qubits. This is highly counterintuitive in the clas¬ 
sical setting, where completely-mixed states cannot be 
used to yield additional information from a phase mea¬ 
surement, and the maximum value for F would be 1 (as 
attained by a single pure qubit). Therefore the enhance¬ 
ment of Fq by l is an entirely quantum-mechanical phe¬ 
nomenon. 

It is tempting to say that the resource that enables this 
enhancement in Fq is the entanglement or discord in er. 
However, a closer look at a in the limit e —> 0 reveals 
that it is an equal mixture of products of eigenstates of 
a x (for which |—) occurs even number of times), 

( a ®* +1 + ° T ' +1 ) - ( 16 ) 


and it is therefore fully-classically correlated [24] , 
Though er is separable, and therefore preparable via unre¬ 
stricted LOCC, it cannot be prepared using bulk LOCC 
operation. Without the Cnot gate used in the state 
preparation, which is controlled on a quantum superpo¬ 
sition, the register of maximally-mixed qubits cannot be 
exploited. 


At this point we can ask whether there is any discord 
present in the final state of the circuit. In Ref. 35] it 
was shown that there is no discord in the output state 
of a DQC1 circuit when the controlled-unitary operation 
is Hermitian, i.e. U = W in Fig. [l] (see also Refs. pB] 
and for further details). The unitary operator Uu( r ), hr 
Eq. (14), is Hermitian if and only if w(r) = 0, i.e. when 
(j) is known to perfect precision. Therefore it may be ob¬ 
served that the circuit in Fig. [2] contains discord for all 
runs except when </> is fully known. Repeating this analy¬ 
sis for arbitrary values of l, to, n > 0 shows that the final 
state is always separable, but has finite discord except 
when u>(r) = 0. The only exception is when l = m = 0, 
in which case the final state has no correlations. We may 
conclude that noisy input states lead to discordant out¬ 
put states in our model, which sheds new light on the 
constant level of discord at the output of DQC1 found in 

Ref. [ID]. 

Discussion .— We have constructed a model of quan¬ 
tum metrology, inspired by DQC1, that uses highly- 
mixed states as its enabling resource. Our most surpris¬ 
ing result arises when the register is taken to be fully 
mixed. In this case, the probe state is classical correlated, 
and yet it can only be prepared via a coherent quantum 
interaction due to the bulk operation constraint of our 
model. Whilst there is no entanglement or discord in the 
probe state, we have found that the state at the output 
always has discord, except in the limit of infinite preci¬ 
sion for the phase parameter being estimated. Our model 
then surpasses the performance of a classical setup when 
only one qubit in the register has a finite amount of pu¬ 
rity. In this case the probe state also has entanglement, 
which is widely understood to be essential for achieving 
precision beyond the SQL. 

Our results provide support for both entanglement 
and discord as enabling quantum resources in quantum 
metrology. Perhaps more importantly, our model shows 
how a large ensemble of highly-mixed quantum systems 
can be of great utility for quantum sensing. Since our 
model only requires bulk coherent operations on the en¬ 
semble, it has the potential to enable a scalable quantum 
technology could challenge state-of-the-art classical sen¬ 
sors in the near future. The biggest practical weakness of 
our model lies in the fact that if even a single qubit is lost 
between the first and last controlled gates, all sensitivity 
is lost — a problem which is shared by any measurement 
device using pure GHZ states or NOON states in the 
context of interferometry [12l . 
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Appendix A: State evolution throughout the protocol 

We begin with the initial state 


_ 1 /|O)(O| 0n <g> pf m (E> If' |O)(O| 0 " O pf m (8) 1 f'\ 
a ° ~ 2 V|O)(O| 0n (8> pf m (8> If' |O)(O| 0 " (8) pf m (8) if7 ' 

After the first Cnot we have 

1 f |O)(O| 0ri (8> pf m <8> l 0 ' |O}(1| 0 " <8) {ptfy x )® m 8 a®' \ 

ai ~ 2 ; +! V|l)(O| 0n 8 (<J xP eT m 8 of 1 \D(ir ® (< T x p,a x )® m <8 1®7 ' 


(Al) 


(A2) 


Next the phase is encoded, but note that u$ commutes with |0)(0|, p € , and a x p e a x 

1 / |O)(O| 0n ®pf m ® U 0 ' |0)(1 |® n ® {utp e a x u\)® m ® { U4> (j x u\)® 1 

CT2 “ 2'+! [e^ n+1 ^ |l)(O| 0n 8 {u^a x p e u\)® m 8 |l)(l| 0n 8 (<7 x p e <r x )® ro 8 l 0 ' 

Next we apply the second Cnot gate 


T.3 


1 / pf m 8 l 0 ' 

2 1+1 \e l ^ n+1 ^(a x u <f> a x p e u] t> )®' m 8 (cr x u^,a x u^)® 1 


e * (n+1)0 (w0p e cr a: uJ ) (7 a; ) 

pf m 8 


0m 8 {u^,a x U^O x 




Finally we apply controlled-^( r ) as well as vgt r \ on the control qubit. Note that v$i r ) also commutes with 
and a x p e <j x . Measuring the control qubit in the basis of a x gives us 


q ± (r) = (±|ct 4 |±) (l ± Re{tr[e 


-i(n+l)(c/>-9(r)) 


|0)(0|' jn 8 (peU,f > a x u l (p a x vl (r) )^ m 8 (n^u^n^) 0 1 /2 l }}\ (A3) 


„t 


(l ±Re{e < n+v >^ e WHr[p e u, i £7 x M^f7 a; t;J (r . ) ] m tr[M ? if7 x tt^f7 a; nJ (j , ) ]'/2'} 

^ 1 A Re ^ < 




cos(o;(r)) + iesin(o;(r))] cos L (cu(r 


(A4) 

(A5) 



